IceRay: An IceCube-centered Radio-Cerenkov GZK Neutrino Detector 
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Abstract 

We discuss design considerations and simulation results for IceRay, a proposed large-scale ultra-high energy (UHE) 
neutrino detector at the South Pole. The array is designed to detect the coherent Askaryan radio emission from UHE 
neutrino interactions in the ice, with the goal of detecting the cosmogenic neutrino flux with reasonable event rates. 
Operating in coincidence with the IceCube neutrino detector would allow complete calorimetry of a subset of the events. 
We also report on the status of a testbed IceRay station which incorporates both ANITA and IceCube technology and 
will provide year-round monitoring of the radio environment at the South Pole. 

Key words: neutrino detection, Askaryan effect, radio frequency 
PACS: 14.60Lm, 95.55.Vj, 98.70.Sa, 84.40.-x 
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1. Introduction 

Continued progress in the determination of the ultra- 
high energy cosmic ray (UHECR) spectrum above 10 17 
eV has established the presence of the Greisen-Zatsepin- 
Kuzmin (GZK) suppression [l[ , resulting from the interac- 
tion of UHECRs with the cosmic microwave background 
(see Fig. [J). Such interactions lead to a "guaranteed" flux 
of UHE neutrinos, although the characteristics of the flux 
depend on the details of the source distribution, UHECR 
composition, and other currently unknown factors. Mea- 
surement of the GZK neutrino flux would not only shed 
light on these issues, but also could indicate the UHECR 
sources themselves via the direction of the individual neu- 
trinos. 

An array to detect UHE neutrinos via their coherent ra- 
dio emission in a dense medium was originally described 
by Gusev and Zhelezhykh 0, based on theoretical work 
by Askaryan. Since then, significant experimental work 
by the RICE collaboration Q has established many of the 
fundamental characteristics of radio transmission in the 
polar ice, and the ANITA balloon experiment [1] has cur- 
rently set the best limits on UHE neutrino fluxes. Direct 
observation of coherent radio emission using an ice target 
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Figure 1: World ultra-high energy cosmic ray and predicted cos- 
mogenic neutrino spectrum as of early 2007, including data from 
the Yakutsk [jj, Haverah Park Q, the Fly's Eye 0, AGASA @], 
HiRes [lj, and Auger collaborations. Data points represent dif- 
ferential flux dI(E)/dE, multiplied by E 2 . Error bars are statistical 
only. GZK neutrino models are from Protheroe k Johnson |lO| and 
Kalashev et al. 
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at SLAC has also confirmed the theoretical foundations 
described by Askaryan [121 ]. 

The expected flux of GZK neutrinos is nevertheless quite 
small, requiring a large array (or in the case of ANITA, a 
huge target) to ensure reasonable event rates. We propose 
to extend the IceCube neutrino detector [l3| to energies 
from 10 17 to 10 20 eV with a sparse array of radio anten- 
nas (IceRay). An initial array of 50 km 2 is designed to 
provide event rates of 0(few/yr) and establish the base- 
line flux level, while a final target array of 300 to 1000 
km 2 could provide O(100) events per year. Centering the 
array around IceCube allows a subset of the events to be 
detected in coincidence, providing complete calorimetry 
of both the initial interaction (via the radio emission) and 
the outgoing lepton (via the optical emission). While rare, 
such events provide a valuable means of cross-calibration 
and reduction of systematics in the absolute energy scale. 

2. Design Considerations 

We discuss several of the design considerations for a 
large-scale radio array in the ice, in particular the operat- 
ing frequency and geometry. 

2.1. Operating Frequency 

We are initially bounded in operating frequency to the 
region between several MHz, where backgrounds may be 
prohibitively large, to around 1 GHz, where the ice be- 
comes opaque. The coherent emission from an Askaryan 
pulse has a peak field strength which rises linearly with fre- 
quency, but the received voltage at an antenna is inversely 
proportional to frequency, so the direct dependence on fre- 
quency cancels when considering the signal-to-noise ratio 
(SNR). However, a dependence on the bandwidth remains; 
specifically, we find 



SNR (X -Eshower 1 



1 GAf 
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for a shower of energy i? s hower, using a receiver with gain 
G and noise temperature T sys , and where Z$ is a reference 
impedance. Therefore, high bandwidth is important, but 
there is no direct dependence on the center frequency of 
the band. 

Other considerations, however, indicate a preference for 
lower frequencies. First, while the peak field strength of 
the Cerenkov emission rises with frequency, the angular 
width of the Cerenkov cone gets narrower |l4| . Effectively, 
this reduces the total solid angle available for detection at 
high frequencies. 

Furthermore, the frequency dependence of the attenua- 
tion length of the ice itself plays an important role. Over 
the 200-700 MHz range, the attenuation length decreases 
by approximately 25-30% [HI]. Because the effective vol- 
ume, to first order, varies as £^ tten , this implies a strong 
loss at high frequencies. The overall conclusion is that a 
high bandwidth, low frequency approach is optimal. Given 



that a bandwidth factor of 5 is reasonably achievable, we 
set a preliminary target frequency range of 60-300 MHz. 

2. 2. Geometry 

Because the radio field attenuation length in ice is of 
0(1 km) one can cover a relatively large area some- 
what sparsely. While deploying detectors on the surface is 
the most cost effective, refraction effects greatly penalize 
the volumetric acceptance. The index of refraction varies 
from 1.79 in the deep ice (below about 200m) to 1.33 in the 
packed snow at the surface [§]. The low-density region is 
known as the firn, and upward-going rays moving through 
this region are bent away from the surface. This creates a 
horizon angle — that is, shallower rays cannot reach the 
detector. This angle gets much less severe as one moves 
deeper into the ice (see Fig. suggesting that deploying 
antennas in holes, say, 50m or 200m below the surface is 
much more efficient. 
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Figure 2: Example of refraction effects for shallower antenna loca- 
tions. Both 50 m (upper) and 200 m (lower) deep antenna locations 
are shown. On the left are the wide-scale ray geometries, showing 
the terminal horizon angle in each case, and on the right the details 
of the ray bending in the near zone are shown. 



IceCube has already developed drilling technology that 
can be utilized for IceRay. While the enhanced hot water 
drill (EHWD) used for drilling the 2.5 km deep holes for 
IceCube string deployment is not mobile enough for our 
purposes, the independent firn drill (IFD) which drills the 
"pilot" holes for the EHWD can be easily moved. The 
IFD is a "hotpoint" -style drill which melts into the firn 
using a cone of closed-loop copper tubing, heated with a 
propylene glycol/water mixture. The IFD currently can 
drill at about 4 m/hour, with an average power usage of 
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approximately 100 kW. The IFD is effective to depths of 
40-50 m, after which pooling water causes power usage to 
spike. Adding a pump to extract this water is a simple 
modification which could alleviate this issue. Ultimately, 
we expect that drilling to 200 m is logistically manageable 
and cost effective, either with a modified IFD or other 
technology. 

2.3. Baseline Configurations 

Given the above design considerations, we focus on two 
geometries for the initial 50 km 2 phase of the array: a shal- 
lower, denser array deployed at a depth of 50 m and with 
36 stations; and a deeper, more sparse array deployed at 
a depth of 200m and with 18 stations. The configurations 
are chosen to have approximately the same cost and vol- 
umetric acceptance in the peak energy region of the GZK 
neutrino flux, around 10 18 eV. Figure [3] shows the sta- 
tion arrangement in more detail. Each station consists of 
three holes separated by 5-10 m, with four antennas (two 
of each polarization, horizontal and vertical) in each hole. 
Directionality is achieved for even single-station events via 
timing information from these local baselines. 
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Figure 4: Volumetric acceptance, in km 3 steradian, of various array 
configurations (d = station depth; TV = number of stations; and 
station spacing), including results from three independent simulation 
chains. 



3. Simulated Event Rates 

The primary IceRay simulation chain is based on Monte 
Carlo code developed for ANITA and SalSA [ll], but inde- 
pendent crosschecks have been performed with ARIANNA 
[171] and RICE simulation chains. The volumetric accep- 
tance of different array configurations is shown in Fig. [4j 
and we note reasonable agreement in the important energy 
range of 10 18 eV. In general, the 18-station deep configu- 
ration gives higher acceptance than the 36-station shallow 
configuration at the higher energies, but drops off at low 
energies due to the increased station spacing. 

Table [T] shows integrated event rates for the two base- 
lines configuration studied. "Standard" fluxes, such as 
ESS with a ACDM cosmology [HI and Protheroe [? ], 
result in approximately 3-10 events per year, but as men- 
tioned earlier, this could vary almost an order of mag- 
nitude in either direction depending upon UHECR com- 
position and source evolution. Iron UHECR models in 
particular tend to produce significantly lower rates fl9( | . 
although these are currently disfavored by measurements 
of the spectral endpoint [lj]. An important point is that 
no irreducible backgrounds are expected, so detection of 
even a few events would be significant. 

One significant motivation to build IceRay at the South 
Pole is to allow for the possibility of coincident, or "hybrid" 
events with the IceCube detector. A or v T event can 
produce both an initial shower and a long-ranged charged 
lepton with the potential for detection in both radio and 
optical channels. A typical geometry for such a hybrid 
event is shown in Fig. [5l 

Such a hybrid event allows cross-calibration of the en- 
ergy scale of either detector, and while such events are 



Table 1: Event rates per year for several classes of UHE cosmogenic 
neutrino models. The "36-50" rates are for the 36-station, 50m-deep 
configuration, and the "18-200" rates are for the 18-station, 200m- 
deep configuration. 



Cosmogenic neutrino model 


36-50 


18-200 




ev/yr 


ev/yr 


Fc UHECR, std. evolution 


0.50 


0.60 


Fe UHECR strong src. evol. 


1.6 


1.8 


ESS 2001,fi m = 0.3, Ov = 0.7 


3.5 


4.4 


Waxman-Bahcall-based GZK-u flux 


4.2 


4.8 


Protheroe and other std. models 


4.2-7.8 


5.5-9.1 


Strong-source evolution (ESS, others) 


12-21 


13.8-28 


Maximal, saturate all bounds 


24-40 


32-47 



rare, they are background-free. Event rates per ten years 
for various GZK flux models are shown in tabled] Adding 
a high energy "guard ring" of strings to IceCube (the 
"IceCube- plus" configuration; see Ref. [13]) increases the 
hybrid event rate by up to a factor of two. We have also 
conservatively assumed here that each detector triggers in- 
dependently; adding sub-threshold cross-triggering would 
also increase the event rate. 

4. Testbed Station 

Recent data from the ANITA flights have demonstrated 
that the South Pole is not a particularly radio-quiet envi- 
ronment (at least in the austral summer), but no capability 
currently exists for year-round monitoring of the natural 
and anthropogenic backgrounds. To understand and char- 
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Figure 3: Left: Baseline 36-station, 50 m depth array, in a plan view (top) and side view (bottom) showing the simulated interaction region 
around the detector. Right: Alternative 200 m depth, 18-station array. 
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Figure 5: Example of a hybrid event where the vertex is seen by four 
surface radio detectors and the resulting lepton passes near enough 
to IceCube to make a detection. 



acterize these backgrounds, as well as to test prototype 
hardware, we have built the IceRay testbed station. 

When deployed, the testbed will be a single surface sta- 
tion with pairs of antennas buried in shallow boreholes 



Table 2: Hybrid event rates for the baseline IceCube, and IceCube- 
plus (1.5 km guard ring), per 10 years of operation, for several classes 
of UHE cosmogenic neutrino models, assuming the IceRay-36, 50m- 
deep radio array. 



Cosmogenic neutrino model 


IceCube 


IceCube+ 




10 yrs 


10 yrs 


ESS 200m m = 0.3, Q A = 0.7 


3.2 


6.4 


Waxman-Bahcall-based GZK-f flux 


3.8 


7.6 


Protheroe and other standard models 


3.8-7.1 


5.0-8.2 


Strong-source evolution (ESS, others) 


10-19 


13-25 


Maximal fluxes, saturate all bounds 


22-36 


30-44 



(2.5 m deep) in the snow. The boreholes are arranged in a 
circle of radius 5 m. Each hole contains a discone antenna 
optimized for vertically polarized signals, and a batwing 
antenna for horizontally polarized signals. Figure [5] shows 
the layout of the station as well as the antenna geometry. 

The data acquisition system (DAQ) combines hard- 
ware elements of ANITA, IceCube, and the digital radio 
modules of AURA (Askaryan Under-ice Radio Array; see 
Ref. [22]). Four antennas (two of each polarization) are 
first fed into a low-noise amplifier chain (with a total gain 
of ~ 76 dB) in a shielded housing. The combined sys- 
tem has a bandpass of 115 MHz- 1.2 GHz. High- and 
low-frequency components are split and separately digi- 
tized with LABRADOR3 ASICs [HI at 2 GSa/s and 1 
GSa/s, respectively, as part of the IceCube Radio Readout 
board (ICRR). The digitized waveforms are buffered and 
transferred through an intermediary board, the TRACR, 
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Figure 6: Layout of the IceRay testbed station. 
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Figure 7: Schematic representation of the data acquisition system 
for the IceRay testbed station. 



to a standard IceCube digital optical module mainboard 
(DOM-MB), which also provides event time-stamping via 
its own digitizer, the AT WD. The DOM-MB communi- 
cates via the standard IceCube communications protocol, 
so the station can be connected to the IceCube cabling 
network and controlled from the IceCube counting house. 
Furthermore, standard IceCube time calibration proce- 
dures can synchronize the DAQ to a GPS clock. Figure [7] 
shows a schematic of the DAQ components. 



range, in order to lower our energy threshold below 10 17 
eV. This would both increase the total event rate and pro- 
vide enhanced opportunities for hybrid events. Further- 
more, other techniques such as acoustic detection of UHE 
neutrinos are developing rapidly, suggesting that a hybrid 
radio-optical-acoustic array may have significant benefits 
for systematics and cost |23[. 

Installation of the IceRay testbed in the austral sum- 
mer of 2009 will allow precise characterization of the noise 
environment and will facilitate further development of the 
50 km 2 array. IceCube construction will complete in 2011, 
and we hope to phase in construction of IceRay at that 
time, as the ability to use IceCube as the core of a GZK 
neutrino detector is an unparalleled opportunity. 
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5. Outlook 

We continue to study open issues with the design of the 
full array (of order 1000 km 2 ), such as power distribution 
and communications. We are also actively pursuing efforts 
to increase the sensitivity of the system down to the cos- 
mogenic kT noise floor of about -114 dBm/MHz, as well 
as working with South Pole station management to con- 
trol the anthropogenic noise in the 60-1000 MHz frequency 
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